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Abstract. Wg introduce a reversible Markovian coagulation-fragmentation 
process on the set of partitions of {!,..., L} into disjoint intervals. Each 

\Q , interval can either split or merge with one of its two neighbors. The invariant 

measure can be seen as the Gibbs measure for a homogeneous pinning model 
|9]. Depending on a parameter A, the typical configuration can be either 
dominated by a single big interval (delocalized phase), or be composed of 
many intervals of order 1 (localized phase), or the interval length can have a 

CU ' power law distribution (critical regime). In the three cases, the time required 

to approach equilibrium (in total variation) scales very differently with L. In 
the localized phase, when the initial condition is a single interval of size L, 
J^ , the equilibration mechanism is due to the propagation of two "fragmentation 

fronts" which start from the two boundaries and proceed by power-law jumps. 
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1. Introduction 

Coagulation-fragmentation phenomena are often modeled by Markov processes 
where the configuration at a given time is a set of "fragments" , each characterized 
by a positive number, to be interpreted for instance as a length or a mass. The 
dynamics then consists in fragmentation events (a fragment breaks into, say, two 
fragments with a conservation of the total length) and coagulation events (say two 
fragments coalesce and the lengths add up). Most of the mathematical literature 
(cf. for instance [2] for a recent review) focuses on mean field models, where the 
rate of coagulation of two fragments is a function only of their lengths ii , £2 and 
does not depend on the two fragments being close or far away in some geometric 
sense. Exceptions are for instance the one-dimensional models studied in [6l[8j[l]: 
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there, fragments are seen as intervals of the real line and only neighboring ones can 
coalesce (in these models, however, fragmentation is not allowed). 

In the present work, we consider a (Markovian, continuous-time) one-dimensional 
process where both fragmentation and coagulation occurs. We introduce the model 
in a discrete setting: on the interval {0,1,..., L} each site contains cither 1 or 
particles (at sites and L a particle is frozen for all times) and fragments arc 
the intervals between two successive particles. Coagulation of two neighboring 
fragments (resp. fragmentation) is then interpreted as the disappearance (resp. 
the creation) of a particle. The transition rates are such that the invariant and 
reversible measure tt^ is proportional to A", with A a positive parameter and n the 
number of fragments, times nr=i(^*)~''~^' where ii is the length of the i*'' fragment 
and p is an exponent in (0, 1). It is well known [5] that there exists a critical value 
Ac such that, as L diverges: i) if A < Ac the typical configuration for tt^ contains 
a fragment of length L — 0(1) and all the others are finite ii) if A > Ac the typical 
fragment is of length 0(1) and the maximal one is of length O(logi) and finally 
iii) for A = Ac the fragments have a power-law tail of exponent p. Here, we study 
how the mixing time (i.e. the time to get close to equilibrium in total variation) 
depends on L in the different phases. In particular, we find that the equilibrium 
phase transition reficcts in a dynamical one: the mixing time is essentially of order 
LP for A > Ac and much smaller (poly logarithmic in L) for A < Ac. 

That equilibrium phase transitions often have a dynamical counterpart is a well 
known fact. This is the case for instance for the d{> 2)-dimensional Ising model with 
free boundary conditions, whose mixing time is of order logi, like for independent 
spins, at high temperature and exponential in L below the critical temperature 
[T^ . In that case the reason for the dynamical slowdown is that below the critical 
temperature the phase space breaks into different valleys (pure phases) separated by 
high energy barriers. In contrast, for our model at A > Ac there is no multiple- valley 
phenomenon or coexistence of phases. Instead, the reason for the Lf behavior is that 
in the A > Ac phase, starting with the very far-from-equilibrium configuration with 
just one fragment, equilibration occurs via the propagation of two "fragmentation 
fronts" which move from the boundaries to the bulk of the system and proceed by 
power-law jumps (and therefore faster than ballistically) . See the discussion of this 
heuristics in Section [3. II 

The reason why this model shows a non-trivial phenomenology despite its one- 
dimensional nature can be understood as follows. It is of course possible to interpret 
the model as a one-dimensional non-conservative interacting particle system (or 
spin system). However, in our case the creation-destruction rate at x depends 
not only on the particle configuration in a finite neighborhood of x but on the 
location of the first particle to the left and to the right of x. In other words, jump 
rates are very non-local. A related fact is that, in the localized phase, equilibrium 
particle-occupancy correlation functions decay exponentially (cf. Theorem [1]) but 
the equilibrium in the center of a finite box A C {1, . . . , L}, conditioned on the 
configuration in {1, . . . , L} \ A, can depend very strongly on the conditioning even 
if the box is large (cf. Rcmark[T]below) . If this were not the case, the "fragmentation 
front" heuristics would fail and the mixing time would be of order log L even in 
the localized phase like for usual one-dimensional, finite-range, non-conservative 
particle dynamics. 
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In the recent works [4] [5] was considered the heat bath dynamics for a (1 + 1)- 
dimensional polymer model, i.e. the path of a one-dimensional simple random 
walk, interacting with a defect line whose invariant measure, when projected on 
the set of polymer-line contacts, is also of the form (|2.1[) with p = 1/2. In this 
context, the A > Ac (resp. A < Ac) phase is called the localized (resp. delocalized) 
phase, a terminology that we will use in the rest of the paper even if we forget 
the underlying polymer model. Despite this superficial similarity, however, the 
two dynamical problems are very different: in [H [S] updates correspond to local 
moves of the polymer configuration, while (in the polymer language) the updates 
we consider in this work are very non-local and consist in performing the dynamics 
directly on the set of contacts. In particular, for the polymer local dynamics the 
mixing time in the delocalized phase turns out to be polynomial in L and larger 
than in the localized phase, due to a subtle metastability phenomenon [H |5] . 

Finally, let us mention that in [IF the question of the dependence of the mixing 
time on the system size was analyzed for another coagulation-fragmentation process 
(the "discrete uniform coagulation-fragmentation process" ) which is however of 
mean-field type (but results there are more refined than ours). 

In the rest of the paper, we abandon the coagulation-fragmentation language 
since it is more practical to think in terms of particles and holes. 

2. The model 

2.1. The equilibrium measure. Let K{-) be the probability measure on the 
positive integers defined by 

K{j)=CKr^'^''\ J = 1,2,... 

with < jO < 1 and Ck a positive constant such that X]i>i ^(j) — 1- 
For L e N and A > let tt := tt^ be the probability measure on 

^L = {V = (7?o,m, ■■■,VL)e {0, l}{"-i'-'i} : ,yo == r/i = 1} 



defined by 



Xniv) "(") 



^^^'^^^zziAyH^^^^+i-^j) (2.1) 

where the configuration rj is identified with the set {xq = < xi < . . . < x^u^^^i = 
L} of sites X occupied by a particle (i.e. r]x = 1). The number niji) is the number of 
particles strictly located between and L in the configuration 77. For Ac {0, . . . , i} 
and 77 6 Hl, we denote by 77A G {0, 1}"^ the trace of 77 on A. For any probability 
measure /i on fi/,, the marginal of fj. on {0, 1}^ is denoted by p,\\. 

It is possible to express the partition function Zl{X) = X^jien ''^iiv) i^i ^ more 
compact way as follows. Let S = {5*0, Si, . . . } be a renewal process on the integers 
with Sq = and inter-arrival law K{-), i.e. P{Si — Si-i = j) = K{j); denote P its 
law and E the corresponding expectation. Then, one has 

Zl{X) = E [A^"-"'(^^^>l{Les}] , (2.2) 

where l{jgs} is the indicator function of the event {3i > : Si = j}. One therefore 
recognizes in Zl(X) the partition function of a homogeneous pinning model with 
pinning parameter log A [HI Chap. 2]. 
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The system undergoes a phase transition at A = Ac := 1: for A > 1 (locahzed 
phase) the partition function grows exponentially with L and the set of sites where 
7^2; = 1 behaves like a renewal sequence with exponential inter-arrival law; for 
A < 1 (dclocalizcd phase) the partition function tends to zero and the set of sites 
where i]x = 1 behaves like a transient renewal process (in particular, the number 
of particles is a geometric random variable). In the next theorem we collect a few 
well-known equilibrium facts [9l Chap. 2]. 

Theorem 1. 
• Localized phase. For A > 1 one has Zl{\) ^ exp{LF{X)) with F{X) > 



IS 



(the free energy). Moreover, the law of {0 < x < L : rj^ = 1} under -k 
the same as the law of {0 < x < L : x £ S} under P\{-\L £ S), where S is 
a renewal sequence of law F\ on the integers with 5*0 = and inter-arrival 
law 

Xa(j):=Pa(5.-^.-i=j)= ^-^^ 



Em e-^(^)™A^(m) ■ 

Correlations under n^ decay exponentially fast in space: there exists c(A) > 
such that for every 0<a<b<c<L 

l^LiVb = l\Va = Vd = l)- nlim = 1)1 < e-^(^)"""((''-'^)'(''-'')) (2.3) 

and 

7ri(r/x = for every a < x < b) < p-'^Wib-^) , (2.4) 

• DelocaHzed phase. For A < 1 one has Zi{X) ^ c{X)K{L) and the law 

of {0 < X < L : rjx = 1} under tt^ is the same as the law of {0 < x < L : 
X £ S} under Fx{-\L G 5"), where S is a renewal sequence of law F\ on the 
integers with So ~ and inter-arrival law 

A'a(j) := Px{S^ - 5,_i = j) = XK{j). 

The renewal S is transient, i.e. P\{Si = cxd) = 1 — A > 0. Moreover, there 
exists c(A) > such that for every 1 < x < L — I, one has 

7ri(77. = l)<c(A)(^-^^A_^ ' (2.5) 

and 

T^l{n{ri) >k)< c(A)exp(-A:/c(A)). (2.6) 

• Critical point For A = 1 one has 

TT 

Here and in the following, c(A) denotes some positive constant depending on A 
which is not the same at each occurrence. Also, considering S as the set of the 
renewal times we write L G S for the event {3i > I : Si ^ L}. 



Remark 1. Based on (j2.3[) one could be tempted to think that in the localized 
phase the marginal at site x of tt^ conditioned to the configuration rj outside a box 
{x — £,..., X -|- ^} would depend weakly in the conditioning, for i large. This is 
false. For instance, it is easy to see that, if (- ■^ log L and we condition on r/y ~ 
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for every y ^ {L/2 — i,L/2 + £}, the probability that 77^/2 = 1 tends to zero with 
L — > CO, while the unconditioned 7r^(77L/2) "is bounded away from zero. 

2.2. The dynamics. The integral of a function / witli respect to a probability 
measure fi is denoted by fi{f), the covariance between the functions / and g by 
fi{f ; g). The continuous-time dynamics {r]{t)}t>o we consider (heat bath dynamics 
or Gibbs sampler) makes updates which consist in adding or deleting one particle 
at the time. Its generator is 

L-l 

Cf^Y.^Q.f-f], /:Ol^M, 

x=l 

where {Qxf){i]) — "^iif \ 'fe' V t^ ^')- I^ words, the dynamics is described as follows. 
Each site 1, . . . , i — 1 is equipped with a clock whose rings form i.i.d. Poisson Point 
Processes of intensity 1. If the clock labeled x rings at some time t, we replace 
r]x (t) by a new value sampled from the equilibrium distribution conditioned on the 
instantaneous configuration outside x at time t. It is a standard fact that tt^ is 
reversible for C. 

If we denote by rj^ the configuration 77 where site x has been flipped, 

(viz = (1 - lx=z}Vz + lx=3(l - Vx), (2.7) 

the generator can also be rewritten in the more explicit form 

L-l 

i£f){v) = E [(1 - '?-)^-('') + ^^-dAv)] [fivl - fiv)] (2.8) 

x=l 

where the "creation rate" c^irj) and the "destruction rate" dx{vi) are given by 

Cx{v) = 777 ^ , if Xi<x<Xi+i, (2.9) 

A + 



K{xt+i - x)K{x ~ Xi) 



and 



dx{ri) = -— r, if x^x,. (2.10) 

K{xi+i-Xi)K{x^- Xi-i) 

K{xi+i ~ Xi-i) 

We recall that in the previous formulas Xi,i = 0, . . . , N{r]) + 1, denote the positions 
of the particles in the configuration 77. More generally, if A is a subset of {1, . . . , L — 
1}. the heat bath dynamics with only sites in A updated and boundary condition 
r G flL outside A is defined by the generator 

^X/ = E[2./-/], /:«X^K, (2.11) 

x£A 

where i^\ ^ {rj E Hl '■ rji^o = ta^}. Here and in the following A'^ stands for 
the complementary set of A in {1, . . . , L — 1}. Of course the probability measure 
^l(' 1 1^" ~ '''A'^) is reversible for C\. 
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2.3. Relzixation and mixing time. The semigroup generated by the Markov 
process {r]{t)}t>o with generator C, see (|2.8p . is denoted by e*^ and its Dirichlet 
form is given by 

x = l tj^Ul 

The corresponding spectral gap is defined by 

gap= inf %^^ (2.12) 

/:O.^M 7ri(/;/) ^ ' 

where the infimum is taken over non-constant functions. The relaxation time Tiei 
is defined as the inverse of the gap and, for any / ; Ql -^ K and t > 0, we have 

Hence Trci measures the speed of convergence to equilibrium in L^(7r^)-norm. 
Another natural and widely used way to measure this convergence is with respect 
to the total variation distance. If /i, v are two probability measures on a finite 
probability space E the total variation distance between ^ and v is defined by 

\\^' - '^ll ■= ^ E l'^(^) - ''(^)l = ™P l^(^) - ^(^)l = inf ^(^ ^ ^)' (2-13) 

c 
where the infimum is taken over all couplings of /x, v with X ^ ^ meaning that the 

law of the random variable X is fi. The mixing time T„iix is defined by 
T,„ix = inf{t > : sup ||/i^ - tt^H < (2e)"i}. 

where fi^ stands for the law at time t of the process starting from the initial config- 
uration a. The choice of the numerical factor (2e)~^ is irrelevant (any other value 
smaller than 1/2 would be essentially equivalent) but with this definition we have 

sup ||/.^-^^||<e-L*/T».xJ. (2.14) 

a-eUL 

With respect to the spectral gap, the mixing time is much more sensitive to 
atypical initial configurations (w.r.t. equilibrium). We have the following general 
bounds between Tj-oi and T„iix 

Trcl<T^ix<log|^JTrcl (2.15) 

A/ 



where n* = min^gn^ t^l^v) 



3. Results 



Our first theorem gives the correct order of the mixing time (up to logarithmic 
corrections) in the localized phase A > 1. 

Theorem 2. Let A > 1. There exist positive constants Ci(A), C2(A) and C such 
that for L>2 

Ci{\)LP < T„,ix < C2iX) LP [log Lf. 
The constant C depends only on the exponent p in the definition of K{-). 
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That the cquihbrium phase transition at A = 1 has a dynamical counterpart is 
imphed by our second result, which shows that in the delocalized phase the mixing 
time grows (poly)-logarithmically in L. 

Theorem 3. Let e > 0. There exists C3 < oo and, for every A < 1, positive 
constants Lo(A, e), Ci(\) such that for every L > Lo{X, e), 

(1 - e) log L < T,„ix < (log Lf' (3.1) 

and Tfci < C'4(A) for every L. The constant C3 depends only on the exponent p in 
the definition of K{-). 

In the critical case A = 1 we can only prove non-optimal bounds both for the 
relaxation time and for the mixing time, which however are sufficient to show that 
the relaxation time has a different scaling than in the delocalized phase: one has 

C5 LP < T„,ix < CeL^+P (3.2) 

and 

C5 LP < T,ei < CeL'+P (3.3) 

for some positive constants C5 , Cq (cf . Section ^ . 

A very interesting open problem is to understand whether Tioi diverges with L 
in the localized phase. 

3.1. A bit of heuristics and comments. The reason for the different scaling 
of the mixing time according to whether A ^ Ac can be understood as follows. 
Take A > Ac and start the dynamics from the configuration without particles in 
{1, . . . , L — 1} (i.e. there is just one fragment), which should be so to speak as 
far away as possible from equilibrium. Then, particles (i.e. fragments) start to be 
created at the two endpoints of the system and the two "fragmentation fronts" tend 
to invade the whole interval {1, . . . , L}. The reason for the front phenomenon is 
that, even if A is large (which favors the presence of particles), creating a particle 
inside a very large set of empty sites is very unlikely (see (|2.9p ) and therefore 
particles tend to be created close to positions where there is a particle already. We 
will see that the position of, say, the front moving rightwards can be approximated 
by a continuous time, positive-jump random walk whose transition rate from x to 
X + d is approximately l/d^~^P. It is then clear from classical results for sums of 
heavy-tailed random variables that the two fronts will meet in a time of order LP. 
Once the two fronts meet, the system is essentially at equilibrium. Making this 
heuristics rigorous is the core of Section [5l One technical difficulty is that the front 
position is not necessarily increasing (due to coagulation events near the front), so 
the front itself will be defined in a suitably coarse-grained sense. This is the reason 
for the appearance of the logarithmic factor (which we believe to be spurious) in 
Theorem [1 

The situation is very different in the delocalized phase. Again start as far away 
from equilibrium as possible, which however in this case means starting from the 
configuration with particles at each site (L fragments). Then, particles start dis- 
appearing essentially independently everywhere (and not starting from the bound- 
aries) and reappearing less often. This is clear for very small A since the ratio 
creation rate/destruction rate is of order A, see p.9|) - p.l0|) . Moreover, the latter 
ratio decreases very quickly during the equilibrium relaxation because the creation 
rates become smaller and smaller as new empty sites appear. The mixing time. 
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then, should not be very different from the first time when all the Poisson clocks 
associated to the sites in {1, . . . , L} have rung, which of course is logarithmic in L. 
In reality, however, proving that mixing occurs quickly is much harder. A natural 
idea would be to apply simple-minded path coupling |101 Sec. 14.2], trying to prove 
that in a time of order 1 the dynamics contracts the Hamming distance between 
configurations. This works well for A very small, as the reader can check, but has 
no chance to work up to A = 1. The next idea would be to replace the single- flip 
dynamics with a block dynamics where one updates the particle configurations in 
blocks whose size depends on A and becomes large as A /* 1. In this case, path 
coupling works but the problem is then to compare the mixing time of the single- flip 
dynamics with that of the block dynamics. For non-conservative attractive particle 
systems with finite-range fiip rates, it would be easy (applying the so-called Peres- 
Winkler censoring inequality, see Section[5]) to show that the ratio of the two mixing 
times depends only on A and not on L. In our case, however, due to the non-local 
nature of the fiip rates, this does not work. Therefore, we had to devise a different 
iterative strategy to prove fast mixing in the delocalizcd phase, that we think can 
be of independent interest. The drawback is that we get a sub-optimal upper bound 
on the mixing time (polylogarithmic instead of O(logi)). As a side remark, com- 
paring spectral gaps for single-flip and block dynamics is instead rather standard 
and this is the reason why we get the optimal result Trci = 0(1) in Theorem [3] 

4. MONOTONICITY 

In this paper a central role is played by monotonicity properties of the dynamics. 
On i7i there is a natural partial ordering: we say that "q < £, '\i r^x '^ S,x for every 
x = 0, . . . , i. Analogously, for A C {0, . . . , L} we write ?7a < Ca if ^x < £,x for every 
x G A. The maximal (filled) configuration ?yx = 1,0 < x < L, is denoted by "-I-" 
and the minimal (empty) configuration rjx ~ Q,x ^ 1, . . . ,{L — 1), by "— ". 

Let A C {0, . . . , L} and r, r' be two boundary conditions. For ^ S Vt\ and t > Q 
let Tj^'^'^it) denote the configuration at time t of the dynamics in A, which evolves 
with boundary condition r, started from ^. Let also /ij be the law of rf'''^{t). 

Lemma 1. Let A C {0, . . . , L}, let t,t' G il^ and ^ G ^A'^' ^ ^A • Assume that 
T < t' and £, < ^' ■ Then, it is possible to couple the dynamics {v^'^{t)}t>o o.nd 
[f]^ ''^ (i)}t>o *^ such a way that, almost surely, rf'''^{t) < 7]^ '"^ (t) for every t > 0. 

Proof. Let Q!<a<a-<6</3. One observes that 

A A 

> 



K(b -a) - K(B - a) 



Kib-x)Kix-a) K{f3 ~ x)K{x - a) 

and therefore 



< 



K{b-x)K{x-a)~ KiP - x)K{x ~ a) 

^ + ^ mf^) ^ + -^ K{P ~ a) 

(With s = a — a,t ^ X — a,u = b — x,v = [3 — b these inequalities are equivalent to 

1 111 

■+ < 



S +t U + V t u 
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which is trivially satisfied.) Recalling the definition (|2.9|) for the rate of creation 
of a particle at cc, one sees that above inequalities imply that the rate of creation 
(resp. of destruction) of a particle at an empty (resp. occupied) site x is increasing 
(rcsp. decreasing) with respect to rjy^y ^ x. The claim then easily follows from a 
standard coupling argument. D 

An immediate consequence is that, under the hypotheses of the Lemma, for every 
t one has /Xj ^ nf '^ where, given two probability measures ^,1^ on $7^, we say 
that fj, ^ 1/ {fi is stochastically dominated by i/) if //(/) < v{f) for every function / 
which is increasing w.r.t. the above specified partial ordering of fi^. Letting t -^ co 
we also see that if r < r' 

^Li-\VA^ = ^aO ^ ^i(-|^Ac = r;.). (4.1) 

It is actually a standard fact that, given A, one can construct all the processes 
{il^'^{t)}t>o foi' different initial conditions ^ and boundary conditions t on the same 
probability space, in such a way that rj^'''^{t) < 77^ ''" (t) for every t > 0, whenever 
T < t',£, < £,' (global monotone coupling). In the following, we use this fact 
implicitly whenever we say "by monotonicity...". The law of the global monotone 
coupling will be denoted generically as P. 

5. Localized phase 

Wc first give a rough upper bound for the mixing time Tijiix(A,.T) of the heat 
bath dynamics in an interval Ac{l,...,-L — l}of length £ when a particle at some 
site X e A is kept alive during the time evolution. Hence we consider the heat bath 
dynamics with updates in A\{x} and boundary condition r such that Tx = 1. The 
generator is given by jC-\\f^i, see ()2.11|) . 

Lemma 2. Let a ^ 3 + p. There exists a positive constant C := C(A) independent 
of £ :~ |A| and L such that 

T„,ix(A,a;)<Cr. 

The proof (see Section 15. 6p is based on a geometric technique introduced in 
pn [7] for bounding the spectral gap of reversible Markov chains, cf. also [13 Ch. 
3]. Such technique was applied for instance in |llj to prove that the mixing time 
of the Glauber dynamics at inverse temperature /? for the two-dimensional Ising 
model in a 771 X n rectangular box is upper bounded by m'^^ exp(c2(/3)n) if n < m. 
The additional difficulty in our model is that the interactions are not finite-range 
and decay slowly with distance (in contrast with the Ising model case) and that the 
transition rates are not bounded away from zero. 

We start with the proof of the mixing time upper bound in Theorem [2] and 
postpone the lower bound to Section 15.31 By monotonicity it is easy to show (cf. 
for instance the proof of Eq. (2.10) in [13]) that 

sup ||/.^-^i||<2Lmax(||Ai+-^^||,||Air-^Lll), (5-1) 

so we get the desired bound if we show that max [\\fJ,'^ — 7r^||, ||/x^ ^ ^lII) ^ {4:eL)~^ 
for some t ^ 0{LP{logLf). 

Since we are in the localized phase we expect that equilibration occurs faster 
when starting from the full configuration "-I-" than from the empty one "— " . Indeed, 
the next result says that starting from "-f " the variation distance from equilibrium 
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is smaller than l/(4eL) after time of order (logi)^+" with a the exponent in Lemma 
[2] (just choose £ ~ ^1/(1+") in the statement of Lemma|3]). 

Lemma 3. There exist positive constants C := C(A),c := c(A) depending only on 
A such that, for any £ < L, 

+ -7r;^|| <Ci(e-=^ + e-=*/^°). (5.2) 



\\N -t^lI 

Proof. This is rather standard, see for instance the proof of Corollary 1.9 in [13] . 
Assume for simplicity that £ is an even integer. If {il'^''{t)}t>o denotes the evolution 
started from the equilibrium distribution, we have by definition (|2.13p of variation 
distance and then by monotonicity (with P the law of the global coupling) : 

L-l 

'{t)) < Y. [nvtit) = 1) - ^iiv. = 1)] • (5.3) 

x=l 



llA^^-^ill 


< nv^it)^ 


For £/2<x<L- 


- £/2 we have 


nvtit) = 1) - 


-7ri(7?, = 1)< 



''+(%/2(i) = 1) - nnrU/2 = 1)] + e-^^ (5.4) 

where P^'+ denotes the law of the dynamics in {1, ...,£— 1} with boundary condi- 
tion rjQ = rji = 1, started from the maximal configuration. We used the exponential 
decay of equilibrium correlations (cf. p.Sp i to deduce \Tr'^{i]x = 1) ^ '"'£ (%/2 = 
1)1 < exp(— c^) and monotonicity of the dynamics to deduce F{ri^{t) = 1) < 
P^'+ (77^/2 (i) = 1) (recall that x is at distance at least £/2 from both endpoints 
of the system). By Lemma [2] (with £ = L) the mixing time T,jiix(£) of the dynamics 
defined by ?*"'+ is 0(£") and we recall (cf. ^A^) that 

\P''Hm/2it) = 1) - Trfi^j,/, = 1)1 < e-L*/Tn,.WJ. 

The terms in (|5.3p where x is within distance £/2 from or L are treated similarly. 

D 

The real job is to estimate the time needed to reach equilibrium starting from 
the minimal configuration. We turn to this problem now. 

5.1. The censored dynamics. Wc start with some definitions. Let 

£:=£{L) = l{\ogL)^+^\, 

/3 being a positive fixed parameter and assume for simplicity that L = (r + 1)£ for 
some r > 2. 

Definition 1. For j ~ 0, ...,?' + 1, we define pj ~ ji and we note Aj = {pj — 
£, . . . ,yj + £} n {1, . . . , L}. For later convenience, we define also "reduced boxes" 
A' as follows. For j = 1, . . . , r let 

A'j := {x e kj : (1/10)£ <\x~yj\< (9/10)^} 

while for j = and j ~ r + 1 let 

Ay.= {xeAr. \x-y,\<{9/l0)£}. 

We partition the box {0, . . . , L} into the r overlapping boxes Aj, j — 1, . . . , r. 
Note that yj is the leftmost point of the box Aj+i. We also introduce a sequence 
of deterministic times T^ = fc(l + £"+^), fc > 0, with a the exponent in Lemma [5] 

Let n'f be the law at time t of the censored dynamics {??'(i)}t>o obtained from 
the following scheme. 
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• The initial condition is the minimal configuration "--" . 

• For t G [Tk, Tfe + 1), the dynamics {f?t}t is given by the heat bath dynamics 
on {1,...,L- 1}. 

• For t G [Tfc + l,Tk+i) the dynamics {ry^jt is obtained by the heat-bath 
dynamics by keeping only the updates in {1, . . . , i — ^}\{yjT J — li ■ • ■ > ''}■ 

Observe that tt^ is invariant for the censored dynamics. 
Lemma 4. We have for every t > 



llA^t 



< 



1/4 -^ri I 



and ^'f. :< /i( ; moreover the functions jjl^ /tt^, fJ-t/T^L '^^^ decreasing. 

This is an immediate consequence of the so-called "censoring inequalities" by 
Peres and Winkler which say that, for a monotone heat bath dynamics started 
from the maximal or minimal configuration, censoring (i.e. eliminating) certain 
updates (according to a schedule which does not depend on the realization of the 
process) increases the variation distance from equilibrium. 

Theorem 4. (HJ, see also [HI Th. 2.5]) Let n eN, < to < ti < ...tn = T and 
Ai C {1, . . . , i — 1}, i = 1, . . . , n. Let fiQ be a law on Q^ such that the function 
rj I—)- iif){rj) / 'K^{if) is increasing (resp. decreasing). Let ht he the law at time T 
of the continuous-time, heat-hath dynamics in {!,..., L}, started from fiQ at time 
zero. Also, let vt he the law at time T of the modified dynamics which again starts 
from fiQ at time zero, and which is obtained from the above continuous time, heat- 
bath dynamics by keeping only the updates in A^ in the time interval [ti_i,ti) for 
i = 1, . . . ,n. Then, 



Imt-tt^II < \Wt 



(5.5) 



and fiT d: i^T (resp. vt ^ I^t)', moreover, the functions [Lt It^\ and vt It^\ are both 
increasing (resp. decreasing). 

Another very useful consequence of monotonicity is the following lemma, which 
ensures that, as soon as equilibrium is reached in some box, the system remains 
at equilibrium forever there. The total variation distance ||/i — i/||a in a box A C 
{0, . . . , L} between two probability measures [i, v on $7^ is defined as the variation 
distance between the corresponding marginals on A; |l/i — i^||a := ||m|a — ^|a||- 

Lemma 5. Let A C {0, . . . ,i}. The functions t M- \\ij.f — tt^Ha, t h^ 



\f^t 



-^^,11 A 



and t I— > ll/ij — tt^IIa are monotone non-increasing. 

Proof. Let us give the proof for instance for ||/i( — tt^JJa. For lightness of notation 
we write here n :~ n^. For 77 G ^II write 77 = {rjAjTjA"), so that for instance 
T^iVA) = t^\a(j]a)- From Theorem |4] we know that the function ry 1-^ ^((7/)/7r(7/) is 
decreasing for every given t. We deduce that the same holds for the ratio of the 
marginals of the two measures: 



M't(a) / M^a) 



< 



Indeed, 



fJ.'ti^A.lW) 



whenever ^^ < ^ 



A- 



(5.6) 



< TT 



7I"(Ca,??aO 

l^tiCA,VA'') 



VA = £,A 



. T^i^A^VA-) 



VA = 6 



< TT 



MK^A'^aO 



. T^iCA^VA- 



VA =^'a 
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In the first inequality we used the fact that /ij(-,77A':)/7r(-,7?A=) is decreasing and in 
the second the fact that TT{g{'riA'')\'riA = •) is decreasing if .g(-) is decreasing (this can 
be seen as a simple consequence of monotonicity of the dynamics in A'^ with respect 
to the boundary conditions in A). From Theorem |4] it easily follows that, if s < t, 
then /i', :< /ij and as a consequence also /i^|a ^ fJ'tlA- The claim of Lemma [5] then 
follows from the following result, applied with fi — /i', |a, i^ ~ ^J■t\A and p ~ ttJa- 

Lemma 6. (^ [14|, Lemma 16.7]J Let p, /i, ly be probability measures on a finite, 
partially ordered probability space. Assume that p ^ u and that the function 
rj ^ ^{if)/ pij]) is decreasing. Then, Wv — p\\ < \\p — p\\. 

O 

The usefulness of censoring in our case is shown by next result which says that, 
if J}y = 1 at some time Tk + 1, then just before time Tk+i the system is close to 
equilibrium in the domain A' . cf. Definition [T] 

Proposition 1. Fix some < j < r + 1 and let lo G IIl be a configuration such 
that Uy. ~ 1 . There exist positive constants C, c (independent of j and £) such that 

\\p'd-Win + l) = uj)- TrilU^. < CI (e-^(*-^'=-i)/^° + e-^') (5.7) 

for every t e [Tfe + 1, Tk+i] . 

Note that, by definition of {Tk}k and t = 1{L), for t = Tk+i the r.h.s. of dSJ]) 
is 0(exp(-c^/2)) = 0{L-P) for every p > 0. 

Proof. We have 
\\p',i- I ri'iTk + l)=uj)- Trill^, < \\p[{-WiTk + 1) = lo) - tt^- \ ry,. = w,.Vz)||^, 

3 3 

A „A/ 



Fl -^^Li■\Vy^ ^^y^'^i)\ 



A',- 



From the definition of A', the fact that ujy^ = 1 and the exponential decay of 
equilibrium correlations (cf. Theorem[T]) we have that the second term in the r.h.s. 
is 0(exp(-cf)). 

The process {ri'{t)}t in the time interval [Tk + 1, Tfc+i) is the heat bath dynamics 
with sites {yi}i not updated. In particular, conditionally on the event {r]'{Tk + l) = 
w} with ujy^ = 1, the particle on site yj is alive during the time interval considered. 
We denote also {'fj~^{t)}t (resp. {'fj~{t)}t) the heat bath dynamics with updates 
only in Aj\{yj}, which starts at time T^ + 1 from the configuration which is "+" 
everywhere in {0, . . . , L} (resp. from the configuration which is empty everywhere 
in {1, ... ,i — 1} except at yj). Let {?7'^(i)}t>Tfe+i be the evolution which starts 
at time T^ + 1 from the distribution 7r^(-|?7yi = ujy-Mi) and where sites {yi}i are 
not updated (note that the law of rj'^{t) is 7r^(-|?/yi ~ ^yyi) for all t > Tk + 1). 
By monotonicity, \\p^{-\T]'{Tk + I) = ui) — 7r^(-|?7yi = Wj,;Vz)||a'. is upper bounded 
by the probability that ri'{t) and 77" (i) do not coincide in A'- (conditionally on 
v'iTk + 1) = cj). Also, it is obvious (under the global monotone coupling, again) 
that {ri'{t))A'. 7^ (^'^(*))a'. implies (?7^(0)a' ¥" ('7~(i))A'. • Then, thanks to a union 
bound over the sites of A' , one deduces 

\\p't{-W{Tk + l)=uj)- n^^{-\rjy^ = ^y.Vz)||A;. < P((ry+(t))A^ ^ (?~r(i))A:) 

< Yl iA^fe)-Arfe)i 

xeA'. 
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where p.^ denotes the law of fj"^ (i) . 

Note that the equihbrium measure of 77+(i) (resp. of ?7^(t)) is n^{- \ 77A';u{ai} — 1) 

(resp. 7r^(- 1 ??{i,...,l-i}\a = ^^Vvj ~ 1)) ^-iid by Lemma [2] the mixing time of both 
dynamics is 0(£"). Therefore, 



|/i+fe)-A-fe)| < Ce-^(*-^^-i 



)/r 



+ hLi-\r]A-u{y,} = l)-7i"i(-|»7{i,...,L-i}\A, = 0, ?7w = 1)IIa^. 



Thanks to the exponential decay in (|2.3p . the last term is 0(cxp(— c^)) (here wc use 
that the points in A^ are at distance at least £/10 from A' ) and the claim follows. 

D 

5.2. The auxiliary process. We introduce now an auxiliary process {C{t)}t>o 
where ({t) = {Ci{t), . . . , C,-{t)) £ {0, l}^^--'^^ defined from the censored dynamics 
as follows. The initial condition is Cj{0) = for every j and (,{■) is constant in the 
intervals [Tk,Tk+i). The value of Cj(^fc) is 1 if and only if 

• either Q{Tk-i) = 1 

• or Cj{Tk-i) = and 77' (T/c-i + 1) = 1. This implies, via Proposition [TJ 
that the dynamics rj'it) is almost at equilibrium at time Tk in A'-. 

We also set by convention Co(0 = Cr+i{i) = 1 for every t > 0. Observe that 
{C{t)}t>o is a non-Markovian dynamics, but the joint process {(C(^)7'?'(*))}t>o is. 
Note that in the process {C{t)}t particles never disappear. The next result says 
that, when there is a particle at j for (^{t), with high probability there is a particle 
for r]'{s) in Aj, for every s between t and a much larger time, say L. 



Lemma 7. There exists c 
every k such that Tk < L 



c(A) > such that for every j G {0, 



1} and 



P [Qin) = 1 and 3s e [Tfc,i], 7^'{s)\^^ = O] < (1/c) e-^'. 

Proof. Let Ah be the event Ah = {r]'y^{Th-i + 1) = 1} and r = mi{h < k : 
Ah occurs}. We have the following disjoint partition 

{Cj{Tk) = l} = Uh<k{T = h}. 

Since A' C Aj and k < L ^ exp(c^/2), the desired claim follows if we show that, 
for every h < k, 



T = h and 3s e [Tk, L], r/(s) |a^. = 
The l.h.s. of (|5.8p is upper bounded by 

3s e [Tk,L] such that 77'(s)|a'. = 
3s £ [Tk,L] such that r/(s)|A/ = 



<(l/c)e- 



-rf 



(5.8) 







^h 



7^'y^{Th-i + l) = l 



where the equality comes from the Markov property of the process {if{t)}t- Let 



H 



1(7/(.s)|a' =0)ds. 



Tfc 
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H >0\ v'yXTh-1 + 1) = 1 < exp{-c£). One has 



(5.9) 



so that we need to show that 
via Markov's inequahty 

H>0\r^'y^ (T,,_i + l) = l]<p[0<if< cM-ai) \ <■ (Th-i + 1) = 1 

+ e'^^E[iJ|4(T,_i + l)-l], 

with a to be chosen. The first term is easily shown to be exponcntiahy smah in £: 
this is an immediate consequence of the Poisson distribution of clock rings, because 
< H < cxp{—ai) requires that there are two clocks that ring within a time interval 
exp{—a£) and recall that exp(— a£) = o{L~p) for every p. On the other hand, the 
second term equals 

p(,/(s)|a^. =0|r;;^(T^_i + l) = l) ds. 

Thanks to Proposition [1] and LemmajSlwe know that, conditionally on 'q'yiTh^i + 
1) = 1, for every t > Tk the system is within variation distance cxp(— ci£) from 
equilibrium in the region A', for some ci := ci(A) > 0. On the other hand, at 
equilibrium one has (cf. (|2.4p ) 

7r(r/'|A;. = 0) < exp(-C2£) 

for some positive C2 := C2(A). In conclusion, the second term in (|5.9p is upper 
bounded by 

and the claim follows choosing a smaller than min(ci,C2). 

D 

It is clear that after some random but finite time (as long as L < oo) one has 
Cj = 1 for every j < r. The following result quantifies such random time. 

Proposition 2. There exist positive constants c :~ c(A) and C := C{p) such that 
for R{L) = [LPf^\ we have 

P(3je {!,..., r}, O(T«(i)) = 0)<(l/c)e-^^. (5.10) 

Proof. Define the stopping time 

r = inf{s > : 3j < (r + 1) such that 7?'(s)|a^ = and Cj(s) = 1} 

and the new Markov process {{C,{t), fj'{t))}t which is constructed like {{C,{t), rf{t))}t 
except that all moves, which would lead to ?7'|a' = for some j such that at the 
same time Q = 1, are discarded. Obviously, the two processes can be perfectly 
coupled until time T- Since from Lemma [7] we see that P(T < L) < exp(— c(A)£), 
it is enough to prove (|5.10p for the new process. 

The process {C(i)}t>o can also be defined by the following formula 

O(Tfe) - O(Tfe-i) + (1 - CATk-i))f}'y^{n-i + 1). 

Let hjik) — max{n < j : C,n{Tk) = 1} denote the position at time T^ of the 
rightmost particle of C, before site j. Then, 

nf}'y,{Tk-i + 1) = l|(C(rfe_i),r)'(rfc_i))) > CO ^ — - (5.11) 
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for some positive constant cq. Indeed, in the whole time interval [Tk-i,Tk-i + 1] 
there is a particle in the box A'^ ,j,_j^> for the process ij' . If ?)(,. {Tk^i) = 1, then the 
l.h.s. of (|5.1ip is lower bounded by the probability that the clock at yj does not 
ring before T^-i + 1- If on the contrary 77' . (Tfc_i) = 0, then a lower bound is given 
by the probability that the clock at yj rings only once and that a particle is created 
there. Given the form (|2.9|) of the creation rate, the fact that K{-) is decreasing 
and given the definition of A'- one obtains the claim (jS.lip . 

Let {C{t)}t>o be the Markovian process, constant on time intervals [Tk,Tk+i) 
where C(t) = (Co(0. Ci W> ■ • • ) ^ {0, l}^"-i'->, constructed as follows: 

• the initial condition is Cj(0) = 0, j > 1 and Co(0) = 1; 

• we set Co{t) = 1 for every t > 0; 

• to obtain C,{Tk) given C,(Tk-i), sample Bernoulli random variables Bj{k) G 
{0, 1}, j = 1, 2 . . . , independent for different j, which take the value 1 with 
probability K{j ~ hj{k — 1)) where 

Co is the constant appearing in (|5.1ip and of course hj{k) = niax{n < j : 
CniTk) = 1}. Then, 

0(r,)=0(r,_i) + (l-0(Tfe_i))B,(fc) .7 = 1,2,.... (5.13) 

Observe that K{-) depends on A and £. Note also that, while {C{t)}t is a particle 
system on {0,...,r}, {C{t)}t is a particle system on {0,1,...}. In both cases, 
particles never disappear once they have appeared. 

The reasoning leading to (|5.1ip shows that {C{t)}t dominates stochastically 
{C(t)}t so that 

P( 3j e {1, . . . ,r}, 4(Tfl(i)) = 0) < P(3j e {1, . . . , r}, CjiTR^L)) = 0). 

In Theorem [5] of next section we will prove that the probability in the r.h.s. is 
0(exp(— c^)), c := c(A) > 0, which concludes the proof. D 

5.3. Front propagation for the process {C(t)}t>o- In this section we study the 
process C(t)- First we give a result about the propagation of its rightmost particle, 
which is at position zero at time f = 0. Then, we estimate the time it takes for 
C(t) to fill entirely the interval {1, . . . ,r} where we recall that r = L/£.{L) — 1 = 
L/L(logL)i+'^J -'1. 

We associate to the process {C{Tk)}k=o,i.... the front position, denoted by /(fc), 
defined by 

f{k) = sup{j > : ^,iTk)^l}. (5.14) 

It is easy to see that {/(fc)}fe>o is a renewal process on {0, 1, ... } with /(O) = 
and inter-arrival law P(/(fc) — f{k — 1) = j) =: Q{j) given by 

00 00 

QU) = K{j) n (1 - k{i)), J > 1, Q(o) = n(i - ^«) (5.15) 

i—j-\-l 2—1 

with K{-) defined in (|5.12p and we recall that it depends on £ = £(L). Observe 
that J2T=o QU) = 1 (just note that Q{j) = A{j + 1) - A{j) if j > and A{j) = 
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Lemma 8. There exists ^o ^ 1 ind S > such that for T{n,i) ~ \nP£^^'^\ we have 
V^(i) > 4, Vn > 2, P(/(T(7i, £)) e (n/4, 3n/4)) > S. (5.16) 

The restriction n > 2 is just to guarantee that there is at least an integer in 
(n/4, 3n/4). The proof boils down to showing the convergence in law of n^^ f{T{n, i)) 
to a non-degenerate stable law of parameter p as n, ^ — > oo. This is rather standard 
and details can be found in Appendix [X] 

Theorem 5. For any e > there exists c{e) (that can depend also on X) such that 
if T{L) = LP£^+P+^ we have for L large 

P(Vj- £{!,..., r}, C,{T{L)) = 1) > 1 - e-=(-)^. 

Proof. Let ^ > ^o a-s in Lemma [51 Then, for some positive S we have 

p(y?e{r/4,...,3r/4}, Q{T{r,£))^o) <{l-S). 

Using monotonicity and Markov's property of the process ^ (monotonicity be- 
ing guaranteed by the fact that K{-) is decreasing) it follows that at time t = 
T{r,i)£^^^, with probability greater than 1 — (1 — 6Y we have (x{t) = 1 for 
some site x € {''/4, . . . , 3r/4}. Conditionally on this event, we can repeat the argu- 
ment on the two intervals {0, . . . ,x} and {x, . . . , r} to claim that we can find two 
sites y S {x/A, . . . , 3x/4} and z G {a; + (r — a;)/4, . . . ,x + 3(r — a;)/4} such that 
Q,{y) = Q,{z) = 1 for t' = T{r,(.)t+^' + T(3r/4, £)^i+^ The probability that aU 
this does occur is bounded below by 1 — (1 — 5Y — 2(1 — SY 

We repeat this procedure n times, until r(3/4)" < 2 (note that n = O(logr) = 
O(logL)). After a time t = E"=i '?'((3/4)-'V,^)^^+^ and with probability greater 
than p=l- YTj^i 2^(1 - ^Y^^' we have Cx{t) = 1 for every x G {0, . . . , r}. The 
proof is concluded when one observes that t = 0(r''£^+''+^) = 0{LPp^'^) and 
p > 1 — U{1 — 5Y > 1 ^ e^^^ for some positive c which depends on e and on 5 
(and therefore on A). 

D 

5.4. Proof of Theorem [2] (upper bound). In view of (|5.ip and Lemma[31 it is 
sufHcient to show 

\\^Jii-^l\\<l/{^eL) (5.17) 

for some t = [^L^ i\og L)'~^ \Tx where we recall (see Section [5?T|) that 

Tl =:0((l0gi)(l + '3'(" + l)). 

By monotonicity, like in (|5.3p . 

L-l 

ii/i*- - ^iii < E ^^l{v. - 1) - nv^.{t) = 1)) • 

x=l 

We consider first the points x in R := U^JIqA'. Thanks to Lemma |4] we have 

E i^^'i- = 1) - ip('?-"(^) = 1)) ^ E W(^- = 1) - ip(^.^(^) = 1)) 

<imax|[/xJ-7r^||A'. 

3 ' 
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From Proposition [2] and the definition (. = [(logL)^+''J we know that, with high 
probabihty, 

C(Li''(logL)(i+^)^J)^l. 
Via Proposition [1] and Lemma [5] this imphes that at time t the dynamics rj'{t) is 
within variation distance exp(— c^) from equihbrium in each box A'-. This shows 
that the first sum in (|5.18|) is o{L~p) for any p > 0. 

The contribution from the points in i?' = {1, . . . , L — 1} \ i? is treated similarly: 
one has simply to redefine the points yj as yj ~ [j + l/2)£ (instead of yj = j£) in 
the definition of the censored dynamics {ri'{t)}t. This way, the union of the new 
boxes U^JIj A' covers R' and the argument proceeds similarly as for R. 

5.5. Proof of Theorem [2] (lower bound). We will prove that for times much 
smaller than L'' the dynamics started from the empty configuration is still far from 
equilibrium (in particular, that there are still very few particles at distance of or- 
der L from the boundaries). The idea is to compare the true heat bath dynamics 
{r]~{i)}t>o to the dynamics {'fj{t)}t>o (which again starts from the empty configu- 
ration) defined as follows: 

• the destruction rates dxir]) in (j2.10p are set to zero; 

• as soon as a particle is created at a site x we put a particle also at the 
position L ~ X. 

Let jit be the law of 7?(i). The standard coupling allows us to realize the two 
processes so that fi{t) > r]{t) for every t. Since \\fj.^ . — 7r^|| < (2e)~^ and 
^L iv\{L/i,...,L/2} = O) ^ with L — )> cxo (we are in the localized phase), we have 
for L sufficiently large 

A^T_(^l{V4,...,V2}^0)>l/2. (5.19) 

Since fj{t) > r]{t) for every t we have /iT„ix('7l{L/4,....L/2} ^ 0) > 1/2. Wc define the 
front position f{t) associated to fj{t) as 

fit) = max{z e {0, . . . ,i/2}, %(t) = 1} = L - ram{i e {L/2, ...,L}, f,,{t) ^ 1} 

where the equality holds since the process f}{t) is by construction symmetric around 
L/2 at all times. Hence for L sufficiently large we have 

P (l- V~(T,„ix) e [1/4, 1/2]) > 1/2. (5.20) 

One can easily prove (see a sketch at the end of this section) that f{t) is stochasti- 
cally dominated by a random variable y{t) such that n~^y{n'') converges for n -^ cxo 
to a one sided p-stable law whose density is denoted g. Let ol = Tmix/i'' and as- 
sume liminfi_^oo a^ = 0. We have 



limsupP (i-\/(T„,i,) G [1/4, 1/2]) < limsupP (ia^^O '/(aL^") > 



L ^ ' L 



4 



< lim sup lim sup 




/•OO 

lim sup / g{x)dx — 

A-S.OO Ja/4 



lA/i 

which contradicts (|5.20p . Hence lim inf a^ > and the mixing time lower bound is 
proven. 
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It remains to give the limit behavior of f{nP)/n. Note that {/(i)}t is a Markov 
process with values in {0, . . . , L/2} and positive increments. The rate of the tran- 
sition from value f to f + x (with x < L/2 — /) is given (cf. (|2.9p ) by 

2 ^ <m. (5 21) 

A + K{L - 2f)/{K{x)K{L - 2/ - x)) " x^+p ' 

where the factor 2 is due to the fact that, if a particle appears at x it appears also 
a.t L — X for the process {fi{t)}t. Therefore, f{t) < y{t) where {y{t)}t is a Markov 
process on {0, 1, ... } with jump rate / -^ f+x given by the r.h.s. of (|5.2ip . Finally, 
the convergence in law of {y(n'')/n}„ to a p-stable variable is immediate. 

D 

5.6. Proof of Lemma [2l If a; is not one of the two endpoints of A, then the 
dynamics in A n {y : y < a;} is independent of the dynamics in An{y : y > x}. We 
can therefore assume without loss of generality that x is the left endpoint of A and 
moreover by translation invariance we let x = 0, so that A\{a:} = {1,...,^ — 1}. 
Let also £ + M = mm{y >£ -.Ty = 1}. 

(2e\ 
Recall that T^ix < log — Tici- Here tt* is the minimum over a G {0, l}{i.-"^ i} 

of 

r{a) ~ ^X\{,}(ct) = ^l (r/|{i,...,f-i} = a I r/y = 0,^ < y < £ + Af, J^,+M = l)(5.22) 
It is easy to see that 

-logTT* <Ci(A) x^, (5.23) 

hence T^ix < C2(A) x tx T,.oi. To prove dESS]), given a e {0, 1}^^— '^^-^l let d 
denote the rightmost particle of a before site ^ (set d = if there is no particle). 
One has then from (j5.22p and from (|2.ip : 

„ ,^, j^(M + £-d)e--'(^)'' 

TT (cr) > c(A) — ^-i 

Tfd~=oK{M + i-d')Zd,{\) 

> ciX)'S^±A^e-^-i^)^ > c(A)i^(^e-^'(>)^ 
- ^ ' K{M + l) - ^ 'k{M + 1) 

since K{-) is decreasing and Zl(X) grows at most exponentially in L, see Theorem 
[H The ratio in K{M + (.)/K{M + 1) is clearly lower bounded by C x (.^^^p and 
this concludes the proof of (j5.23p . 

We use a path combinatorics argument to get an upper bound on the relaxation 
time. For any 77,77' e 9.\ = {r/ G 9.l ■ ?7{i,...,f-i}<= == T{i,. ..,£-!}<=} let x^, i = 1, . . . ,75 
denote the set of ordered sitesin{l,...,^— 1} such that rjx i^ Vj^. Consider the path 
7(7^, ry') in the configuration space, of the form (770, . . . , ?7p-i) with 770 = r/,r]p = 77', 
Vi+i — vT (<^f- (|2.7p ). The number p is called the length of the path and is denoted 
by 17(77,77')!. Plainly, one has 17(77,77')! < L Let A be the set {(77, 77^); 77 G r2£, 1 < 
z <f\ and for e = (a, a^) G yt let 

Q^'^ - ^^^^ "^ [ d'{a) if aUl. (5-24) 
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By Theorem 3.2.1 of [15] we have 

^'■°'-Ti{o7;T ^ \liri,v'WivW{v')}- (5.25) 

7(r;,T;')9e 

Consider for definiteness the case CTz = 1, the other case being essentiahy identical 
by reversibihty. One has that dz{a) is bounded away from zero, uniformly in a, L, t, 
so it remains to prove that 

From the way the path ^{rj^rj') was constructed, we see that ^{rj^rj') 3 (cr, cr^) 
implies that rji coincides with ai for i > z while r][ = ai for i < z and r/^ = = 1— ctz. 
One can identify the configuration a with the positions of its particles up to position 
£ + M, = I'o < • • • < x'^(„),i ~ i + M and similarly for 77,77'. Call a the index 
such that x"^ = z, b the index such that x^' = z and c the largest index such that 
x^l < z (observe that x]! = x'^_i). One then sees that 

7r^(77),r^(77') _ 1 iti^i^^Kix]^^ ^)] IV;^ [XK j x]'^, - xf )] ^^^e) 






< ^^^. m^i^Kix]^i--]m^i.'i^K(x]^,-^ )] 



(5.27) 



where 



n((T) 

z^=Y. n [^^(4+1-4)]- 

aeni j=Q 

Since K{-) is decreasing, the expression in ()5.27p increases if we replace x]! by z. 
Finally, it is immediate to see that 

y^ u'=oi^K{x]^i - ^]mK{x^c'+i - ^)] n;Lilm{xf^i - ^f)] 

is equal to TT'^{r]z = !)<!, which concludes the proof. D 

6. Delocalized phase 

In this section, we prove Theorem[3] The first subsection is devoted to bounding 
the mixing time and the second one to bounding the relaxation time. 

6.1. Mixing time. The lower bound Tmix > (1 — e)logL is pretty easy. Let 
1 <C flL -C i and let S be the first time such that all clocks labeled ql < i < L — a^ 
have rung. Since, by (j2.6p . 7r^(77i = for every a^ < i < L — ai) — > 1 when 

L — > cx), Oi — > 00, one has 

WfJ-i - Trill > 7ri(7?. = V fli < i < L - Oi) - P(77+(t) = V ol < ^ < L - a^) 

> l + o{l)-F{S <t). 

(6.1) 
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Note that S has the same law as the maximum of L — 2aL independent standard 
exponential variables so that P(5 < (1 — e) logi) = (1 — L^^^^^^)^^^°-^ which goes 
to whenever e > and the mixing time lower bound follows from (j6.1|) . 

In order to get the upper bound on Tmix in Theorem [31 we generalize a bit our 
model. We define fig := {0, 1}^ and, given r € ^z, ^ an interval in Z. we consider 
the dynamics in A with boundary conditions r on A'^. If r^. = for every x the 
dynamics is defined as the limit as n — > oo of a dynamics with boundary condition 
r'") which equals 1 for |a;| > n and otherwise (the limit exists by monotonicity of 
the creation/destruction rates w.r.t. the boundary conditions) and the equilibrium 
measure is concentrated on the empty configuration. Recall the notation /i^'^ from 
Section |4] for the law of the dynamics started from a, with boundary condition r. 

The proof of the mixing time upper bound is based on an iterative procedure 
(on the size of A) , together with a path coupling argument by Bubley and Dyer [3] 
(see also [lOl Sec. 14.2]). 

Given two configurations ?/, 77' G Oa = {OjI} 1 '^^ introduce the Hamming dis- 
tance d{ri,rj') = X^ieA 1*^* ~ v'li ^ 1^1 which counts the number of discrepancies. 
Given two random laws /i, v on J^a, define their Kantorovich distance as 

dK{n,v):^ inf ¥.d{X,Y). 

(X,Y): 

Then one has the following theorem, which is a direct application of a result from 

Theorem 6. Let the interval A and the boundary condition t he fixed. If there 
exists 7 > and t > such that 



(6.2) 

(6.3) 
(6.4) 

r 

(actually T^(A) depends only on |A|, so we will write rather T^(|A|)). 
We will prove the following. 

Proposition 3. Fix 1/(1 + p) < C < 1 and < (3 < C(l + p) — 1 sufficiently small 
and let 

e{L) := [c2{X)L'^\ . (6.5) 

There exists Ci^p) and an integer Lq(X) such that for L > Lo{X) one has 

Tp{L)<C4{p)T^{i{L)). (6.6) 

Let us show how this implies the mixing time upper bound 







dxilh 


, Mt ) < e ^ 








for every 


a, a' such that d{a,a') — 1. 


, one has 












dKiPnt ^ l^t 


'")<e-T"d(a,a') 








for every 


n G N and 


every couple (a, 


a'). 








Let us 


fix some /? 


> and define 












Tp{A,T) 


:= inf i t > : 


sup dKiPt""^ 

d{a,a') = l 




")< 


1 
|A|^ 


and 




T^(A) := 


= supT^(A,r) 









A ONE-DIMENSIONAL COAGULATION-FRAGMENTATION PROCESS 21 

Proof of Theorem\^ (upper hound). If l^ \-) denotes the application of the map 
X 1—^ (.{x) n times, one has &'^\L) < Lo{X) for n < c{p)\og\ogL, for some finite 
c{p), if L is large enough. As a consequence, 

T^(i)<5:-(logi)^=('')C6(A) 

where essentially Cq{X) = Tp{Lo{X)). 

Since the maximal Hamming distance between configurations in A is |A| and 
thanks to Theorem[6]one has, for the dynamics in A and irrespective of the boundary 
condition t and on the initial conditions cr, cr', 

if n > 1//3, where /? enters the definition (|6.3I) of Tp{A, r). On the other hand, for 
every t > and a S I^a, denoting for simplicity fi^'^ by /i^ and by n the equilibrium 
measure in A with boundary condition r, one has 



WfJ-t -^W < s^pdK{p-t,fJ-t 



Indeed, 



iia^^-^ii = Je 



< sup ll^if - Mt' 



n' 



77,77' 



By definition of the total variation distance and Markov's inequality, one has 

U1 - M?' II < nx ^Y)^ P(d(X, y) > 1) < E(d(X, Y)) 

for every couple of random variables {X, Y) such that AT ~ /i^ , y ~ /ij' . Choosing 
(X, Y) as those which realize the infimum in the definition of Kantorovich distance 

between pi and pi , one one obtains (|6.8p . From (|6.7p we deduce that Ha'S^s ^""11 — 
l/(2e) for every initial condition a, i.e. the mixing time is upper bounded by 
nS. This proves the upper bound in p.ip (actually, a stronger statement which is 
uniform in the boundary conditions). D 

Proof of Proposition\^ Consider cr, cr' G A (with |A| = L; for definitencss. say A = 
{1, . . . , L}) such that d{a, a') = 1 and assume for definiteness that a contains one 
particle more than a', so that in particular cr > cr'. Fix some boundary condition 
T outside A and, for every cr e f^A, denote for simplicity p^'^ by p^. We have by 
definition of Kantorovich distance and monotonicity 

dK(/i^Mf ) < Ed(rf (t),r;-'(t)) = ^ [p^ (r^,) - pf (tj,)) (6.9) 

xeA 

where P denotes the global monotone coupling (so that r]'^{t) > ri°' {t) for every 
t). For simplicity we split the analysis into three cases. Fix some g such that 
C > g > 13/ p and let —a = maxj^; < : Ty = 1}, 6 = minj^; > L : Ty = 1}. 

(i) First case: £{L)/2 < x < L — £{L)/2. Consider the dynamics on A := {a; — 
£{L)/2, . . . ,x + £{L)/2 — 1} (assume for simplicity that £{L) is even and note 
that A C A) with boundary condition fy = l for y ^ A. Call p^ its law at 
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time t when the initial condition is ^ G fi^ and cah tt its invariant measure 
(which is nothing but a space translation of Tr^f^-i). One has by monotonicity 

M^fe) - fifiv.) < /^"fe) < A+fe) < Hv.) + Ut - n- (6-10) 

Now choose i = C4 T^(^(L)) with C4 = C4(p) > 1. The first term in the r.h.s. 
of (|6.10p is upper bounded by c{X)i{L)^^^P (cf. (|2.5p ) and the second one 
by £{L)-'^''f^ (thanks to ^^, the definition of Tp and Theorem [H). These 
bounds are uniform in the boundary condition r. As a consequence, 



(/.;^(7?,)-^^'(,7,)) <^ (6.11) 



«(L)/2<x<L-£(L)/2 



if C4 is chosen sufficiently large and L is large enough (we used that /3 < 

C(i + P)-1). 

(ii) Case max(i9 - a, 1) < x < e{L)/2 (the case max(L9 - (fe - L), 1) < (L - x) < 
£{L)/2 being analogous). This is not much different from Case (i). Note that 
max(L^ — a, 1) <C i{L)/2 since we assumed g < (■ 

Let this time A C A be the interval of length £{L), whose leftmost border is 
max(— a, x — £{L)/2) and consider the dynamics in A, with boundary condition 
Tj, = 1 for y ^ A. Again one has (|6.10p where last term in the r.h.s. is negligible 
if C4 is chosen large enough. As for the first term, one has from p.Sp 

~( )< ^(^) < ^(^) I ^(^) 

^^''' - min((a + x),ey+P " (a + x^+p I^+p ' 

Summing on max(L^ ~ a.,^) < x < £{L)/2 one obtains a quantity which is 
0{i~P)+0{L~sP) = o{L~^) due to the way /? and g were chosen. In particular, 
for L large enough we have 

(Air(^.)-M"'fe)) <^- (6.12) 

max(L9-a,l)<a:<£(L)/2 

(iii) Case < a; < max(L^ — a, 1) (or by symmetry (L — x) < max(L^ — (6 — L), 1)). 
This situation requires more refined equilibrium estimates. The rough bound 
t^tiVx) — ^^t iVx) < /J'tiVx) is not sufficient any more and one needs to show 
that the difference between the two terms is small. 

We can clearly assume that a < L^. We consider this time two dynamics 
on A = {1, . . . , i{L)}: the first has boundary condition Tx = t^ for x < and 
Tx = I for X > i{L) while the second one has boundary condition Tx — Tx 
for X < and ri^^ = for x > £{L). CaU /^^*^'* (i = 1, 2) their law at time t 
when the initial condition is ^ G fl^ and tt*-*' ,i = 1,2 their invariant measures. 
Then, by monotonicity, 

^J'?ivx) - ^J.t ivx) < tjl'p'^iiix) - ^it'^'~(^lx) 



E 



< 



^^'\vx) ^''H^x)] + ll^i'^'^ - -<^'ll + IImP'-^ - -<^)| 



Once more, the last two terms are smaller than £{L)~^'^^ and they arc negli- 
gible if C4 is chosen large enough. It remains to upper bound the difference 
T^ iVx) ~ "^ {Vx)- One has by translation invariancc 

T^'-^Hvx) = 7r^+£(L)+ife+a == MVy = 0, 1 < y < a). 
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Recalling Theorem [1] this can be rewritten as 

(1) _ Eo<^<. Kxja + y)Px{x ~ y)Pxii - x + 1) 
^ ^'^'^ j:o<y<e+iKx{a + y)Px{i-y + i) 

where with some abuse of notation we let Pa("^) '■— Pa("- G S) and we wrote 
£ = £{L) for lightness of notation. Analogously, one can write 

(2). ,^ _ J2o<y<x -^A(a + y)Pxix - y)'Rx{i - x) 



TT^^^fe) 



PxiSi=(x) + j:o<y<iKxia + y)Iixie~y) 

where Ra('T') := Pa(^* > 0,j > n : Si = j) (i.e. the probability that the 
transient renewal 5* contains no points after n) and altogether 

tt'^'Hvx) - TT^^Hvx) = Yl Kx{a + y)Px{x~y) 

0<y<x 

' Px{i-x + l) Rxji-x) 

Eo<y<^+i Kx{a + y)Px{i - y + I) 1 - A + Eo<y<^ ^xia + y)Tlx{e - y) 

(6.13) 

where we used Pa(5'i = oo) = 1 — A. Since the renewal S is transient, it is well 
known pj App. A] that Px{n) ~ jj^jy^Kxin) for n — >• oo. As a consequence, 
one sees immediately that 

|l-RA(n)| =0(n-''). (6.14) 

It follows that 

^^(^"") = I + o(i-n 

l-^ + Eo<y<eKx{a + y)Ilx{e-y) l-XEU^^V^ 
because the denominator of the term on the left hand side is equal to 

a a-\-l oo 

1-^A'a(2;)- J2 Kx{y)[l-Rx{i-{y~a))]- J^ ^>^^y) 

y=l y=a+l y=a+l+l 

and, decomposing the second sum according to y < a + i/2 and y > a + £/2, 
we see this equals 1 — J2y=i ^a(2/) + 0{i^P). 

Now we need the following more accurate estimate, which is proven in 
Appendix IB] there exists w — w{p) £ (0, p) such that 

A 

To upper bound the first term in the square bracket in (|6.13p . we estimate 
separately numerator and denominator. The numerator is 

where we used the fact that {£ — x) > £/2. 
As for the denominator, it equals 

i/2 e+1 

XY,K{a + y)Px{£-y + l) + X ^ K{a + y)Px{i - y + 1). (6.16) 

y=l y=i/2+l 



PA(n) = j^^-^K{n){l + 0(n-)). (6.15) 
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By using ()6.15p and the decreasing character of A'(-), the first sum can be 
bounded below by 



A2 



(1-A)2 
A2 



y=a+l 



(1 - A)^ 



;K{e)\ J2 K{y) + o{rn\{i + o{£-n) 



^y>a+l 



A2 



(1 - A)^ 



;Kii)\ Y, Kiy)\{l + 0{e-n) 



^y>a+l 



and the second sum by XK{i + a + I) X]«=o ^^iv)- 

Note that X]^o ^^ (2/) ^'^ ^^^ average total number of renewal points in S: a 
simple computation on geometric random variables shows that it equals 1/(1 — 
A). From this fact, it is easy to deduce that J2y=o 'P^iv) = l/(l-A) + 0(^"''). 
Observe also that, since we are assuming a < L^, by taking g sufficiently close 
to l3/p and /3 sufficiently smafi one has K{a + 1 + 1) ^ A'(£)(l + 0{e-P)). 
Putting everything together, the quantity in (j6.16p is lower bounded by 

A/(l - \r R'ii) [1 - A ^ R'iy)] (1 + or-)). 
In conclusion, the l.h.s. of (|6.13p is given by 



E 



0<y<x 



Kx{a + y)Px{x~-y) 



1-AE:=iA%) 



K{£-x + l) 
W) 



K{i) 



-i + or"") 
- 1 + or"') 



<(l-A)-i Y. ^^a(2/)Pa(x-2/) 

<i<y<x 

where we used the fact that A'(-) is decreasing and 'Ylm Ki^) = 1. 

Uniformly in the range 1 < a; < £/2 one has, for a positive constant C{p), 

K{e-x + i) ^,^,.x 

Also, the renewal equation gives 

Y Kx{y)Vx{x-y)=V^{x) 

0<y<x 

so that finally 

c(A,p) [^,1 



7:^'>{Vx)~7r^'>{Vx)< 



(2)/ 



3.17) 



(6.18) 



Sumnfing on a; < max(i^ — a, 1) < i, one gets a quantity which is 0(£ "') 
o{L^^) if (3 is small enough. In particular, for L large enough 



max(L^— a,l) 



Y (/i"(^x) -/if ('?.)) < 



8Lf^' 



(6.19) 
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Putting together Eqs. ()6.1ip . (|6.12p and (|6.19p one obtains that the sum in ()6.9|) 
is smaUer than 1/L^ , which finishes the proof of Proposition [3l 

D 

6.2. Relaxation time. Wc prove here the second part of Theorem [3l i.e. Tici < 
C4(A). The proof is decomposed in two steps. First, in LemmaOwe show that the 
Dirichlet form of the initial dynamics is comparable to the Dirichlet form of a block 
dynamics whenever the size of the blocks remains finite, i.e depending on A but not 
on the size L of the big box. Thus, it is sufficient to prove that the relaxation time 
of the block dynamics does not depend on L. In the second step (Lemma [TO|l we 
show, by using the path coupling argument by Bubley and Dyer [3] , that if the size 
of these blocks is sufficiently large then the relaxation time of the block dynamics 
is independent of L. 

The first step would be standard (cf. e.g. [121 Sec. 3]) for, say, the heat bath 
dynamics of the nearest neighbor Ising model, where transition rates Cxio'ja') are 
uniformly bounded from below and depend on cr, a' only in a finite neighborhood of 
the site x where the Poisson clock rings (finite-range interactions). In our case, both 
uniform positivity of the rates and finite-range dependence fails, cf. the definitions 
(|2.9p - (|2.10p of creation/destruction rates. 



The block dynamics is denoted by {fl{t)}t>o- The size of the blocks equals 2^+1, 
where in this section I depends on A but not on L. For any a;£{l,...,L — l}lct 
Ag{x) = {y : \y ~ x\ < i} f) {1, . . . ,L — 1} and A^(x) the complementary set of 
A£{x) in {0, ... , L}. The generator C of {77(i)}t>ois given by 



L-l L-l 



{cfM - Y. [^-/('?) - /('?)] = E E^-(^'^) [/(^) - /(^)] 



x=l 



n,v 



where {Qxf)iv) = ""l [f\VA'{x))- The rate Cx{r],fj) is zero if the restrictions to 
A^(x) of r/ and fy are distinct. If instead 77a<:(2:) = J7a<=(2:)j let a := a{x,vi) (resp. 
h := b{x,ri)) be the rightmost particle of rj in A^(x) to the left of x (resp. the 
leftmost particle of rj in A1{x) to the right of x). Then, 

c- ('?' V) = ^ „,a.b. . (6-20) 

with W"'''(cr) = A*"' n7'=o -^(-^j+i ~ -^j)' ^0 = a, ■Zfc+i ~ b, and Zj,j ~ 1, . . . , fc the 
ordered points z of Ai{x) for which cr^ = 1. In this formula, it is understood that 
for fc = 0, W^^''{a) = Kib - a). If w e Ql, we use the short notation W;^^''{uj) for 

WS'K^Aeix))- 

The Dirichlet form associated to C is given by 

4(/;/) =7ri [/; -Cf] = ^E E ^iiv)cx{v,v) Uiv) - f{v)f ■ 

x — l rf.f} 

Lemma 9. Fo?' any £ > 1, there exists a constant C :— C{i,X) independent of L 
such that for any function f : Ql -^ R, 

sdf-J) < csuf-J). 
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Proof. In the proof, C denotes a positive constant (depending only on X,p and £, 
but not on L) which is not the same at each occurrence. 

For any configurations rj, f) which coincide outside of Af (x), we denote by ^x{'^^ v) = 
{77 = r]^^\ il^^\ ■ ■ ■ , V^^^ = v} ^ sequence of configurations such that ?y('+i) is ob- 
tained from 77^*'' by adding or deleting a particle at some site Xi € Ae{x). The 
sequence is chosen such that p := p{rj^f]) is minimal. Observe that p is at most 
2£ + 1, independently of 77, 77, x, L. 

Write f{fi) — f{ri) ~ X]?=o(/(^^''^^'') ^ /('7^"'^)) ^^d use Schwarz inequality: 

p-i 



P-l ^A 



T^tiv) c^ivJl) 



2 



Here, c{r]'^^\7]'^^'^^'>) stands for the rate of change from 7;^^' to 7;^^+^' for the original 
dynamics {rit}t- By inverting the different sums, we are left to show that for any 
z € {1,. . . ,L — I}, any ui G il^, 



^R{u,z,x,r],'fj)^ E E -*- 



xGAeiz) (r,,ij) xeAt(z) {i^SD LV / A . ; 

is bounded by a constant independent of w. z and L. Observe that the two sums, be- 
cause of the indicator function, can be carried over a set whose cardinal is bounded 
by a constant (independent of aj,z,L). Therefore, it is sufficient to bound each 
term of the sum. Since {uj,ll!^) G 'Jxiv^ v)-, the restrictions to A^(a;) oi u!,u!^,ri,fj are 
equal and a(x,uj) = aix^rf) ~ a{x^f)) and similarly for b. Thus, 

R{uj,z,x,Tj,ii) = l^(„_„.)g^^(^ -)j— — _-— ^--— (6.21) 

c(uj,uj )Wx' (w)ECTe{o,i}'^f(-) W^:r (o-) 

To fix ideas, we assume that w^ = (the other case can be treated similarly). 
The rate c(a;,w^) corresponds then to the creation of a particle at site z. Since 
(w,w^) G Ixin^'n) '^e have rjz =0, fjz = 1- Let u> a (resp. 7; < &) be the rightmost 
particle of uj to the left of z (resp. the leftmost one to the right of z). We have 
from ((2^ 

XK{v - z)K{z - u) 
^("' "^) ^ XK{v^z)K{z-u) + K{v^u)- ^'•'') 

We distinguish two cases according to the positions of u, v with respect to a, b. 

First case: a < u (the case v < b being very similar). This means that u G Ae{x), 
so that K{z ~ u) > K{2£ + 1). Then, from ()6.22p one realizes immediately that 
c(a;,a;^) is lower bounded by some positive constant. On the other hand, one has 
that either 7;„ = 1 or ?)„ = 1 (or both) since the path "fxiViV) is assumed to be of 
minimal length. Assume for definiteness that the first case occurs. Then, it is easy 
to see that W!^'''{r])/W^'^{uj) < C and then ()6.21|) is upper bounded by a positive 
constant (clearly EaG{o,i}A.(=") WS^H'^) > W'°''(^))- 

Second case: u — a and b — v. This means that the restriction of uj to Ai{x) 
is empty. Since (w,w^) G 7a; (t?,??), iiz ~ 1, and we get that ^"'''(77) < CK{z — 
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a)K{h — z). It follows easily (using the definition of A'(-)) that 



c{uj,uj'')W!i'\uj) 



<c 



K{z-a)K{b-z) 1 
K{b - a) ^ A 



<C' 



so that R{uj, z, X, r/,fj) < C also in this case. 



D 



Lemma 10. There exists £ := ^(A) such that the block dynamics {?7(0}t>o has a 
spectral gap bounded below by a constant C := C(£, A) > independent of L. 



Proof. Wc recall that the Hamming distance d((T, w) between two configurations 
a,uj G Hl is defined by ^((7, w) = J2x=i \'^x—'-^x\- By [3] (see also TheoremlHlabovc), 
it is sufficient to find a generator C acting on functions / : fii x fi^ 9 (?/, rj') ^^ M 
which acts like C if f depends only on rj (or only on ?/) and such that there exists 
a positive constant 5 > satisfying 



for every a; S fl^ and 2 G {1, 



, L — 1}. We write 



(6.23) 



L-l 

El 



:^r(/)-/('^,^)] 



where z^^'" is some coupling between 7i'^(-|cr|A<:(a;)) and 7r^(-|a;|Ac(a;)). Let a € VIl 
and z € {1, . . . , L — 1} such that CTz = 0. We know that 7r^(-|o'^|A'=(a;)) dominates 
stochastically 7r^(-|o'|A<:(2-)) a-nd in this case we require v'^''^ to be a monotone 
coupling. We denote by < a < z (resp. z < b < L) the rightmost site to the left 
of z (resp. the leftmost site to the right of z) which is occupied in the configuration 
a. We have (using d(cr, a^) ~ 1) 



mia,an= E [^r^(rf)-i]+ E K'''^w-i 



a::|a; — zl <£ 



li \x ~ z\ < £ then a\A''{x) = c''^|a<;(2:) and we have v^'"^ (d) = 0. Thus the first sum 
is equal to — \{x G {1, ...,£ — 1} ; |.t — z| < €}| < —£. In the second sum, observe 
that if a; > b+£ (or similarly if x < a~£) then after the update of the box A({x), the 
distance between a and a^ is again 1, so that these terms do not contribute. The 
reason is that the presence of a particle in b "screens" the effect of the discrepancy 
in z: one has 7r^(-|cr|Aj(2;)) 
left of the second sum is 



and 7r^(-|a'^|A'^(x)), so the coupling is diagonal. What is 



E 



E 

xeli(a,z,b) yeKt(x) 



{' 



?fe 



'n\K'j,{x) 



lAJ(x) 



riy 



V\a'^{x) -0-|A|(:r) I 



(6.24) 
where /^ (a, z,b) ~ {x £ {1^ . . . ,L ~ 1} : \x ~ z\> £,a ~ £ < x <b + £}. It remains 
to show there exists a constant £q := £o(A, p) such that if -^ > -^o then this sum is 
less than £/A. We restrict to the sum carried over the sites z -^ £ < x <b + £ since 
the other case can be treated similarly. 

By monotonicity, we can and will assume that a is empty in {l,...,z — 1}. 
This way, the first term in (|6.24p does not change and the second one decreases, so 
altogether the sum increases. 
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For any a < /?, we denote by /i"''' the law of the position s of the leftmost 
particle in Ae{x) C {a, . . . ,/3} conditionally to the fact that there is a particle on 
site a and on site /? and no particles in {a + 1, ... /3 — 1} n A'^{x). Note that there 
could be zero particles in A^(a;) C {a, . . . ,/?}, i.e., /i"''^ is a sub-probability. 

For any x such that z+£ < x < b+£, we denote bx the leftmost site y > {x+£)AL 
such that (jy = 1 (observe that bx = b ii x < b — £ — 1). Then we have 



E {- 



Vy 



yeAe{x) 



'yUjCa;) — O- |Ac(:r) 



Vy 



'7|aj(x) = 0-|AJ(a;) I 



^a;+^ 



^ {^il^''- (s) ~ fxl^'- (s)) Trf,^,^] [J2vy\vu^0ioi-x + i<u<bx 




.0,6,. 



seAf(2:) 



^l/=s 



^(.) 



..25) 



where the last inequality follows because, in the delocalized phase, the average 
number of particles in a box is uniformly upper bounded by a constant, ttj^ ^^ , 

denotes the equilibrium measure tt^ -s+i translated by s. Below wc estimate the 
last line in (|6.25p distinguish various cases according to the values of x, bx etc. Recall 
that /i^^^" is a sub-probability, so in many cases (when the box Ai{x) is sufBciently 
far from z,bx) the sum XiseA (x) l^x^''{^) i^ quite small. In the remaining cases, we 
really have to estimate the variation distance between ^^•''^ and /i°'''=" . 

By using the same notations as in the proof of Proposition [31 we can write 



M?'''(5) 



T.s<.<x+, Kxjs - a)Pxiv - s)KxW ~ v) 
A'a(/3 ~a)+ 'Ex-e<u<v<x+e -^a(w - a)Px{v - u)Kx{P ~ v) 



Case 1: z + £ + l<.T<fe — ^— 1. We have bx = b. By using the renewal equation 
(|6T7)) and Pxin) < CK{n) (cf. (l6T5)) l. we get 



E E ^^l'^ 

x=z+e+l seAe{x) 

b-e-1 

x=z+i+l seAf(a;) 



b-e-1 



(^)^ E E 

x=z+e+i seAi{x) 

K{s - z)K{b - s) 
K{b - z) ■ 



'Ls<v<bKx{s ~ z)Px{v ~ s)Kx{b - v) 
Kx{b-z) 



The latter sum can be written as 

b-l 



J2 [{s-z)-^ + {b-s)-'Y^'m{z,s,b) 



s=z+l 

with Ni{z, s, b) the number of sites x such that s e Ai{x) 



s - z if z + 1 < s < z + 2£, 
iVf (z, s, 6) = i 2£ + 1 if z + 2£ + 1 < s < & - 2£ - 1, 
fe-siffe-2^<s<6-l. 
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By using [(s - z)-^ + {b - s)-^] ^^^ < C [{s - z)-^^+p^ + {b - s)-^^+p^] , we are 
left to estimate 

1 x-^ b — s 



E 77^+ (2^ + 1) E 77— uT^+ E 



s=z4-l ^ ' s=z+2e+l ^ ' s=b~2e ^ ' 

and 

Z+2S. , \ 6-2£-l -I 6-1 -, 

a_s)i+P + (2^+^^ 2^ (5_s)i-Kp+ Z^ (fe-s)P 

s=z-|-l ^ ' s=z-|-2£-|-l ^ ' s=b-2i ^ ' 

which are of order i^~P uniformly in z and b. 

Case 2: sup(& — £, z + £+l) < x < b + £ (note that there are at most 2£ + l possible 
values for x). Let e > sufficiently small and assume to simplify notations that 
e^ > 2 is an integer. In the sequel, the constant C^ depends on e while C does not. 
Case 2. a: z < x — £ < z + el: As in Case 1, we have 

2^ Mx'^(s)<C 2^ j^.f^^ _ . ^C ^ [{b^~s) +{s-z) \ 

sSAf(x) si^At(x) seA£(a;) 

which is finite because X]n>i ^ *^ +oo. Since there are at most el 

possible values of x, the contribution to (|6.24p is upper bounded by Cel. 
Case 2.b: x - (. > z + e(. and b^ - {x + i) > eL Again 

s£Af(a;) seAf(2:) 

seAf(2;) 

It is easy to show that summing over the allowed x this term contribute to 
(lOil) at most CJ^^-P\ 
Case 2.c: x — I > z + el and bx — [x + 1) < el: Here at last we use the upper 
bound ||Aif/== -mS'''" II in (E2S1)- We define F{a, s) = N{a, s)/D{a, s) where 



6^ — a , 

/ sKvCx+t 



X — T / brr — u\ 

D{a,s) = l+ Y. Vx{v-u)Kx{bx-v)\l-y—^\ 

x-i<u<v<x+f. \ ^ / 

Note that Mx''"(''^) = Piz,s), /^°'^"(s) = i^(0,s) so that 



-(i+p) 



Jo 
Observe now (using also Z? > 1) that \daF\ < \daN\ + N\daD\. Using 
Pa('^) < CK{n) and the renewal equation (j6.17p we get easily 

/ bx-A-''"-'' 1 1 

N<c\i- -^ — — -- < a 



bx - a) [bx - sf^+P) - '{bx- sY^+P) 

( bx-s\-''^'' 1 1 

\do.N\<c[\-- -— -<C, 



bx-a (^bx~ay{bx-s)P- '{bx~aY{bx~s)p- 
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The last inequalities follows from a < z, bx — s = bx — s + {s — z) > bx — s+si, 
the fact that the function x i—> x/{x + ei) is increasing and bx — s< {2 + e)i. 
Using similar bounds one obtains 

1 1 i^^f 

|9„D| < Ce-, TTj y ^ < a 



Thus, integrating w.r.t. a and using bx ~ z > 2ei, we obtain 



/ \daF{a,s)\ da < 
Jo 



C{e) 



1 e-P 



{bx - s)P {bx - sf+P 



The sum over s S h.g{x) of the r.h.s. is bounded above by C^i^P- Since we 

have now to take the sum over a set of x's of cardinality at most 2£ + 1, we 

conclude that the contribution of these x's in (|6.24p is bounded by C^i'-^^P^ 

Putting everything together, we have proven that (j6.24p is upper bounded by 

Cei'^^P + Ce£. Taking e sufficiently small and then i sufficiently large, this is smaller 

than i/A, which concludes the proof. 

D 

7. The critical point 
Here A = 1; we write tt^ := tt^^^ and (see Theorem[T]) 

Zl := Zl{X = 1) - Ci X j^, Ci > 0. 

It would be easy to prove Tmix > C's LP as we did in the localized phase: starting 
from the empty configuration, for times much shorter than LP the chance of creating 
a particle between eL and (1 — e)L is small. On the other hand, at equilibrium 

T^drix = for every eL < x < {I - e)L) - C2 x e^P , C2 > 

and the mixing time lower bound follows for e sufficiently small. 

What we do instead is to prove that C^LP < Troi < CqL^^p, which implies 
directly the mixing time bounds via (j2.15|) and the fact that, as proven in Section 
EKcf. dElSD), one has log(l/7r*) = 0{L). 

To get a lower bound for Tid, recall its variational definition (|2.12p and choose 
./('?) = J2x=i "^^ {^^^ number of particles between 1 and L — 1 in the configuration 
77). It is well known (and can be seen by an elementary computation) that TTL{f) ~ 
CiLP and ■KL{f\f) ~ C2L^p for some Ci,C2 > 0. Next we compute the Dirichlet 
form. By reversibility (1 — rix)cx{'>])Tr{r]) = TT{ri^)ri^dx{r]^) so that for every / 

L-l 

£L{f-J)^J2 E ^L{v)M,.=i}dx{ri){f{v)-f{vir- 

x=l riGilh 

In our case {f{i]) — ./(?7^))^ = 1 and the destruction rates are upper bounded by a 
constant, so that 

£L{f;f) < C7TL (Y.^V. = i] = C'^Uf) < c'Lp 

with c, c' some positive constants. Altogether, we have proven that the spectral 
gap of the chain is 0{L^p). 
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Finally, the upper bound Troi < CqL^^'^- Note that the proof of Lemnia[2] (which 
holds for every A) gives directly Troi < cL'^^P; we want to improve on this estimate, 
using the fact that wc arc at the critical point. 

From (|5.25p . (|5.26p . the fact that \'j{i],ri')\ < L one sees that to prove Trci = 
OIL^'^P) it suffices to show there exists a positive constant C such that, for every 
(T G f^L and z = 1, . . . , L — 1 satisfying cr^ = 1, 

where we recall that we identify the configuration a with the positions of its particles 
up to position L: = Xq < ■ ■ ■ < 2;^(o.)+i ~ ^ ^^'^ similarly for rj, rj' (we also let a 
the index such that x'^ ~ z, b the index such that x]^ — z and c the largest index 
such that x^ < z). One has, since a;^' = z, 

E X{K{x']^,-x'j)=Z^^cz^-\ (7.2) 

Next, since as we already observed that x^ — 2;^_i, one has 

E n ^K'+i - ^f ) = P(^ - <-i e ^, ^ n {1, . . . , (z - <_i)} = 0) (7.3) 

j)':l(^=0, J=c 



where like in l\2.2\ we let P denote the law of the renewal S with inter-arrival 
distribution K{-). The probability in (j7.3p equals 

L 

^ K{y)V{L -yeS). 

Using the fact that P{L ^ S) ~ Zl{X = 1) ^ cLP~^ one sees easily that 

L ( L"-^ if (z - < 1 ) < L/2 

E KiyML-yeS)<cJ t^^^,. ^ ^ (7.4) 

Putting together (|7.2p and (|7.4p one sees that the l.h.s. of (|7.ip is upper bounded 
by a constant times L'' as we wanted. D 

Appendix A. Proof of Lemma [H] 

We first prove that n^^f{T{n,£)) converges in law to a non-degenerate stable 
law of parameter p. The only point which requires some care is that the law Q{-) 
depends on the parameter i. This will prove the claim ()5.16p for n > nQ,£ > £q for 
some finite rio,^o- For every positive fi we have 

E (exp {-{fi/n) f{T{n, i))}) = [E (exp(-(Ai/n) 1(1)))^^"''^ = [1 + {^if^M - l)f^"''^ 

(A.l) 
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where £((u) is the Laplace transform of the law Q defined in ()5.15p : 

oo 

-« = l-^Q(j)(l-e-^-'"). 
Let us start from the identity 






Q(i) = coAA'((j + 1)^) - CO 



XK{{j + 1)£) + 1 



K{j) 



n (l-^W)-l 



(A.2) 



The dominant term is the first one: one has 



^X((j + f)£)(l 



'M) 



CiifiP ( f°° l-e 



(A.3) 



£l+p 



J.1+P 



■dx + o(/i) 



as /x — >■ 0. The second term in (jA.2p gives a negligible contribution. Indeed, simply 
note that 



As for the last term, it is easy to see that for £ sufficiently large 

c(A) 



n i^-m)-i 

i=3 + l 



< 



(1+pi^j + l)p 



(A.4) 



(A.5) 



(just use that for (. large one has that K{i) is small for every i, so that (1 — K{i)) « 
eyi^{—K{i)) k, exp(— c(A)/(f^"''''(i + 1)"^"'"''))). As a consequence, one sees that the 
contribution to £(/i) from the last term in (|X2|) is 0(^™''^(i'2'')£-2-2p). Altogether, 



£(/x) -l = -A 



gi+P 



1 + 



min((l-p),p) 



£1+P 



o(/i) 



(A.6) 



for some positive constant A depending only on A, p and cq. Plugging this estimate 
into (jA.ip and recalling that T{n,i) ~ [nP£^^''\ we obtain the desired convergence 
for n,i ^ 00. 

Finally, we prove the claim (j5.16p for £ > £o a-nd n < uq. For this, we bound 
below the probability in (|5.16p with the probability that, among the first T{n,i) 
jumps of the renewal {/(fc)}fe, there is one of length n/2 and all the others have 
length 0. In formulas, 

P(/(r(n,£)) G (n/4,3n/4)) > T{n,e)Q{0f'-''-^^-^Q{n/2), 

where the prefactor T{n, t) comes from the choice of the location of the long jump. 
Recalling the definition oi Q{) and the estimate (|A.5|) one sees that this expression 
is lower bounded by a quantity which depends only on n and therefore by a positive 
constant, since n < no. 
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Appendix B. On the Green's function of a transient renewal 

Here we prove the following statement, which is a bit more general than (|6.15|) : 

Proposition 4. Let K{-) be a probability law on the positive integers such that 
n ^^ K{n) is decreasing and K{n) '^ Ck /n^^'' for some p G (0, 1) and some 
Ck > 0. Let A < 1 and K\{-) :~ XK{-). Then, ifPx is the law of the transient 
renewal S with inter-arrival law K\{-), one has 

PA(7VG^)^^°°-^^(l + 0(iV-)) (B.l) 

(1- A)^ 

for some < w < p. 

Without the explicit estimate on the error term this result is well known, see 
e.g. Th. A.4]. 

Proof. Let P := Pa=i. Start by writing 

N 

and remark that if C = C'(A) is sufficiently large, the terms with k > ClogN are 
negligible and give a contribution which is say o{N^^^'^'') = o{K{N)N^'^"). Next 
we show that 

P{Sk=N)^kK{N){l + 0{kyN''")) (B.3) 

for some positive c and some < w < p, which together with the previous observa- 
tion about the large values of k implies the desired bound (|B.2| . The lower bound 
is easy: 

AT 

P(5fe = ^) > kP{S, ^N,Si> (3/4)iV) > kK{N) ^ P(5fc_i ^ N - x) 
= kK{N){l - P(^fe_i > iV/4)} 



> 



> 



kK{N) |l-(fc-l)p[5i > N/{A{k-l))]\ 



In the first inequality, we used the partition of {Sk = N} into the sets {Sk = 
N} n {Si — Si-i > (3/4)A^}, i = 1, . . . , fc, and the exchangeability of the law of 
{Si — Si-i}i<k- The second one follows from that A'(-) is monotone decreasing. 

Let us now prove the upper bound. For any fixed < (5 < 1, by the same 
argument as above, we have 

P{Sk ^N)^k P{{Sk ^N}n {S, >{N- N')}) 

+ Pi{Sk ^N}n{$i<k : S,^ ^,_i > (A^ - N^)}). 

On one hand one has 

N 

kP[Sk=N,Si>{N-N^)]=k Yl K{x)P[Sk-i=N-x] 

<kK{N -N^) ^kK{N){l + OiN^-'^)) 
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(there is at most one jump longer than N — N^ if N is large). On the other hand, 
P(5fc =N,$i: {S, - 5,_i) > (iV - N'))) 

< kF{Sk =N;{N~ N^) > Si > (5,; - 5,_i), i = 2, . . . , fc) 

= k Y, P{Sk=N:x = Si>{S^-S^^i),i = 2,...,k) 

x=(N/k) (B.6) 

<k Y^ K{x)T'{Sk-i ^N ~x)< kK{N/k)P{Sk-i > iV*) 

x={N/k) 

< k{k - l)K{N/k)V{Si > N^/{k - 1)) < ck^KiN) x 0{N-^p). 

D 
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